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An upscattering spectral formation model for the prompt 
emission of Gamma-Ray Bursts 
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ABSTRACT 

We propose a model for the spectral formation of Gamma Ray Burst (GRB) 
prompt emission, where the phenomenological Band's function is usually applied 
to describe the GRB prompt emission. We suggest that the GRB prompt emis- 
sion is mainly a result of two upscattering processes. The first process is the 
Comptonization of relatively cold soft photons of the star off electrons of a hot 
shell of plasma of temperature T e of the order of 10 9 K (or kT e ~ 100 keV) that 
moves sub-relativistically with the bulk velocity Vb substantially less than the 
speed of light c. In this phase, the Comptonization parameter Y is high and 
the interaction between a blackbody-like soft seed photon population and hot 
electrons leads to formation of a saturated Comptonization spectrum modified 
by the sub-relativistic bulk outflow. The second process is an upscattering of the 
previously Comptonized spectrum by the plasma outflow once it becomes rela- 
tivistic. This process gives rise to the high-energy power-law component above 
the peak in the EF(E)— diagram where F(E) is the energy flux. The latter 
process can be described by a convolution of the Comptonized spectrum with 
a broken-power-law Green function. Possible physical scenarios for this second 
upscattering process are discussed. In the framework of our model, we give an 
interpretation of the Amati relation between the intrinsic spectral peak photon 
energy and radiated energy or luminosity, and we propose a possible explanation 
of the GRB temporal variability. 

Subject headings: radiation mechanisms: general — radiative transfer — gamma- 
ray burst: general — stars: massive 
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Introduction 



Understanding the physical processes which give rise to the observed spectra of the 
prompt emission of Gamma Ray Bursts (GRB) is presently one of the most exciting is- 
sues studied by bo th the theoretical and observational community. The Band function 



( IBand et al. 1 119931 ) up to now widely used to describe their prompt emission is a pure 
phenomenological model. It consists of two low-energy and high-energy powerlaws with 
photon index F\ and I^, respectively, smoothly joined at some energy E^. Among the 
physical models proposed for a possible or igin of the Ba nd function it is worth mentioning 
the optically-thin synchrotron model (e.g. iTavani Ill996[ ). in which the electron population 
of the relativistically expanding shell is accelerated by internal shock collisions, eventually 
producing a supra-thermal powerlaw-like distribution in the tail of the Maxwellian. This 
mixed electron population interacts with possible turbulent magnetic fi elds frozen in the 
plasm a emitting synchrotron photons and forming synchrotron spectra (IRees fc Meszaros 
20051 ). The observed break in the energy distribution is thus naturally explained by the 
transition from the optically thick to opticall y thin emission regim e. However the current 
spectral analysis of both time-integrated /e.g.. ICrider et al. 111997 1) and time-resolved (e.g., 



Crider et al. 111998 I ; iFrontera et al. Il2000t iGhirlanda et al. 1120031 ) prompt GRB spectra has 



revealed some problems related to the synchrotron/SSC emission models. 

Many theo retical efforts hay e also been performed in the thermal (photospheric) inter 



preta tion (e.g., iThompson Ill994t iMeszaros fc Rees 



20111 ). From the observational point of view, iRyde &: Pe'er 



2000 



3eloborodov 



ip 

- 



2010: iToma et al. 



(120091 hereafer RP09), analyzed 



a sample of GRBs using archival BATSE data and found that the time-resolved spectra can 
be fit by a high temperature (10 9 K) blackbody (BB) spectrum plus a powerlaw (PL). Their 
model, albeit phenomenological, actually strongly points in favor of the presence of a pho- 
tospheric emission process at the origin of almost 50% of the total emitted energy. In fact, 
the low-energy threshold of BATSE prevented RP09 to use more detailed thermal models, 
but the good fits provided by application of the BB+PL model could point in favor of the 
presence of a photospheric (Compton-saturated) emission plus a second process giving rise 
to the PL-like hard X-ray emission. As also argued by RP09, one of the advantages in con- 
sidering thermal processes is to reduce the kinetic-to-radiation conve rsion efficiency, which 
is difficult to account by other theoretical models. However recently IZhang fc Yan I (j201lf ) 
demonstrate that the observed high efficiency is a drawback of the internal shock model, but 
is not an issue for the models that invoke dissipation of a Poynting flux. 

Another issue related to the GRB origin is the nature of the progenitor. Long GRBs 



are preferably obse rved at high redshift, with an average z > 2 ( jJakobsson et al. 1 12006 



Fiore et al. 



2007 I ). They are concentrated in small, irregular galaxies and show strong 
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evidence of association with Type Ic supernovae (IKelly et al. 1120081 ). Both theoretical and 
observational works point in favor of a Wolf-Rayet star with mass higher than ~40 M 
(IRaskin et al. 1120081 ) and stripped H -envelope as the most likely progenitor of long GRBs 
(the collapsar model, IWoosley 1 1201 ll ). These progenitor masses are expected to leave, in 
their explosion, a central remnant and it is generally supposed that the GRB engine is 
powered by gravitational energy release of a torus of matter debris accreting onto a black 
hole. In some sense, a temporary microquasar forms in the star after core-collapse, and 
consequently torus accretion may power an expanding (relativistic) jet. However, unlike 
the case of AGN or microquasar jets, the environment surrounding the outflow is not the 
interstellar medium but the star photosphere which may play a significant role in dragging 
the jet kinetic energy. The confinement level of the expanding jet depends on the external 
environment; if the condition #jMj < 2 (where 8j and Mj are the jet opening angle and Mach 
number, respectively) is satisfied, the pressure equilibrium with the surroundings and jet- 
structure is almost maintained, while in the opposit e case strong shocks inside the jet may 
form and the overpressure allows a free expansion ( iGehrels et al. 1 120121 ) . Baryon- loading 
may also be a key ingredient in the determination of the jet bulk velocity. 



In this context, it is worth mentionin g the theoretical wo rk by lLazzati et al. I ( 120091 ). 
which used a numerically relativistic code ( IFryxell et al. I l2000h to evaluate the evolution of 
a jet as it leaves a massive progenitor star after core-collapse and propagate s to the photo- 



spheri c radius, where radiation is eventually released. In their simulations, lLazzati et al. 



( 120091 ) considered a 16 M Wolf-Rayet progenitor star evolved to pre-explosion and a jet 
with initial opening angle 8q = 10° and Lorentz factor 7 = 5 at a distance R = 10 9 cm from 
the center of the star. Although limited to this particular configuration, the authors showed 
that the theoretical light curves are in good agreement with the observed ones. In addition, 
the photospheric temperatures are in the range 100-300 keV, in turn leading to values of the 
rest-frame peak energy E p consistent with observations. However, also progenitors associated 
with much more massive stars (Population III stars) cannot be excluded, at least for farthest 
GRBs. Indeed, supernovae associated with these stars with Helium core M c > 100M Q , have 



been observed also in the closer Universe ( IGal-Yam et al. II2009I ). Such progenitors, with 



initial ma ss in the range 130-260 Mm, have been shown to explode due to pair-instability 



(PI) (e.g., lLanger et al. 1 120071 ; IWoosley et al. 1 120071 ; IWoosley 1 1201 if ) with no compact rem- 



nant. Thus GRBs could be the result of helium burning i n degenerated condit i ons, i n which 
the burning occur s in deflagration regime [see details in IChardonnet et al. I (120101 )]. Also 
in this case, (e.g., Ilmshennik et al. I Il999l ; iBychkov et al. 1 120061 ) have shown that a shock 
runs away from the burning (reaction) zone which leads to the development of an outflow 
(jet-like) structure. In terms of this model, the GRB spectrum would be originated from 
Compton upscattering of soft photons (probably < 1 keV) in the hot sub-relativistic bulk 
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outflow region. The hydrodynamical simulations of IChardonnet et al. I ( 120101 ) show that the 
electron temperature of the hot corona region can be as high as T e ~ 10 9 K. 

The rapid temporal variability on time scale of order AT < 1 s and shorter, observed 
during the GRB prompt phase, implies that the sources are compact with size R < cAT ~ 
10 11 cm. To avoid the problem of high compactness, that would imply high optical depth 
for pair production r 77 ^> 1, a hig h Lorentz factor of the relativistic motion of the emitting 



plasma is the standard scenario (IPiran et al. I Il999h . Alternatively, to decrease r 77 , it is 



needed to consider a much more extended volume over which the photon field is distributed. 

In this Paper we offer a model to explain the GRB prompt spectral formation in the 
context of a photospheric scenario in which the main process is the Comptonization of the 
relatively soft photons of the star by a hot subrelativistic outflow within an area close to the 
photospheric radius (optical depth of 3-5) likely symmetric with respect to the rotational 
axis of the star. 

We st art from the theo r etical and numerical results of the Comptonization problem 
reported in iTitarchuk et al. I ( 119971 . hereafter TMK97) and iFarinelli et al. I ( 12008 I . hereafter 
F08), but for the case of an early subrelativistic bulk outflow phase produced during the 
supernova explosion. This physical scenario i s natural in the case of PI-SNe and it is foreseen 



in the case of an electromagnetic outflow (jLyutikov et al. 1 120031 ). while in the collapsar 
model, the jet formed inside the star is assumed to be initially relativistic. Whether this 
condition can hold, may depend on several parameters, such as the initial jet energy, the 
chemical composition of the star envelope (in particular the presence or not of a H-envelope) 



the core angul ar momentum and the influence of magnetic torques (IGehrels et al. 1 12012 
Woosley llioill ). 



In §2] we present the radiative transfer mathematical formalism, necessary to investigate 
both the spectral formation in the subrelativistic phase, and the formation of the high- 
energy PL component. In §3] we demonstrate that our upscattering model of GRB radiation 
reproduces the Amati relation between the EF(E) peak energy E p and the GRB radiated 
energy E iso or luminosity L gr b. In §H we put arguments for the observed high and slow 
temporal variabilities of GRB emission. In §5] we estimate the resulting luminosity of GRB, 
discuss a formation of the high-energy PL component and our scenario in the light of the 
Fermi results. In 36] we draw our conclusions. 
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2. Upscattering of soft photons of the star in the hot plasma (Compton) cloud 

2.1. Spectral formation in sub-relativistic outflow 

The emergent spectrum produced by Comptonization of soft photons off a thermal elec- 
tron p opulation in the presence of bulk motion has been first formulated by Bland ford fe Payne 



Cowsik fe Lee 



1982 



(|198ll. herea fter BP81) and later investigated by several authors (e.g. 
Colpi I Il988h . Depending on the sign of the divergence of the velocity field, the resulting 
spectrum can b e harder (V • V < 0) or so fter ( V ■ V > 0) with respect to the case of 
static medium (jSunyaev fe Titarchuk I Il980l hereafter ST80). Examples of the first case 
(converging flow) a re reported in TMK97 and F08, w hile the second case (diverging flow) 
was investigated by lLaurent fe Titarchuk I (119991 120071 . hereafter LT07). 



The basic idea of the photon scattering effect in the case of bulk inflow and outflow 
is presented in Fig. 1 of LT07: a photon emitted outwards near inner boundary and then 
scattered at a certain point by an electron moving with velocity Vi, is received by an electron 
having velocity V2 with frequency v-i — V\ [1 + (V\ — V2) ■ n/c] where n is a unit vector along 
the path of the photon at the scattering point. In a diverging flow (Vi — V 2 ) • n/c < and 
photons are successively redshifted, until they are scattered to an observer at infinity. On the 
other hand, in a converging flow (Vi — V2) • n/c > and photons are blueshifted. This pure 
geometrical (Doppler) effect can be analytically described when the electron temperature is 
very low, as the contribution of the thermal velocity vector V t h is negligible compared with 
the dynamical (bulk) velocity vector Vt,. If, on the other hand, the plasma temperature is 
high enough so that Vb and V t h become comparable, then the emergent Comptonization 
spectrum is marginally modified by the presence of the bulk motion. A detailed analytical 
treatment of photon diffusion (Fokker- Planck approximation) in a converging flow is reported 
in TMK97 and F08. Below we will start from the same equations treated by these authors, 
but we will consider the case of a fluid which is subrelativistically expanding (diverging flow), 
namely we consider a velocity field with V • Vb > 0. 



2.2. Formulation of the bulk motion problem. The main equation and its 

solution 

In this section, we provide the mathematical formulation of the problem required to 
compute the emerging spectrum due to Comptonization of soft photons by electrons in a hot 
corona bounced to the top of the star and moving outwardly with subrelativistic velocity (see 
Fig. [1]). The spectrum can be derived using the Fokker-Planck expansion of the radiative 
transfer equation in the presence of bulk motion (BP81, see also Eq. 13 in TMK97): 
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dn 
~dt 

v 2 dv 



+ V b ■ Vn 



-Vn 



kh 



m P c 



n + 



L3A; 

kT e + m e V b 2 /3 dn 



1 dn _ ^ T 
+ 3^V-V b + 



h 



dv 



+ j(r, v), 



where Vb is the bulk velocity field of the outflow, kT e its electron temperature, k(r) = 
N e {r)<JT is the inverse of the scattering mean free path, n = n(r,u) is the zero-moment 
occupation number of the specific intensity and j(r, v) is the source term. 

Let N(r) = Ar (ro/r) M be the radial number density profile of the bulk outflow, with 
radial outward velocity given by 



/3(r) = M^) 2 ~" 



(2) 



which is easily derived from the continuity equation in spherical geometry. Here (3 = Vq/c is 
the dimensionless outflow velocity with respect to speed of light at the bottom of the shell. 
Thomson optical depth of the flow from some radius r to infinity is given by 



r 



N e a T dr' = r (r /r) 



n-i 



(3) 



where To is the optical depth at the bottom of the shell. Using the given radial density profile 
together with definitions (j2J) and (j3J), the steady-state Fokker-Plank (diffusion) equation flTJ 
becomes 



1 ^ ^ <9n ^ 1 <9 

T dx x 2 dx 



x 



i I f dn , 
fbTr- + n 

dx 



3(t, x) 
Qk 



(4) 



where 



5 = — KU ( — ] 
~ 360* -1) r W 



(3-2^/(^-1) 



(5) 



In equation (j4]), L T is the space operator (see iTitarchuk et al. 1 120031 ) . x = hv/(m e c 2 ) 
and = kT e / (m e c 2 ) are dimensionless photon energy and electron temperature, respectively, 
while f h = 1 + (H/c)7(30). 

We note that S in equation (jSJ) is constant only for u=3/2, which corresponds to the 
free-fall velocity profile case (BP81), while for /x ^ 3/2 it is a function of the optical depth 
r. For the particular case of constant outflow velocity (/i = 2), equation (TJJ becomes 



1 





1 d 2 n 



3dr 2 



3r or 



3r 9a; cte 



<9n 

/bo" + " 

ax 



0A; 



(6) 



-7- 



The dependence of 5 on r does not allow to find a solution of the equation ([6]) with the 
method of variable separation (space and energy). However, LT07 demonstrate that for a 
constant velocity profile, 5 can be replaced by its value at some effective optical depth r cff , 
which is a fraction of total optical depth r of the expanding shell. Thus, replacing r with 
r e ff in equation (jSJ) we can now rewrite equation (j3J) in the form 

(C T + L x )n(x,T) = -^-s(r), (7) 

where L x and £ r = _1 L T are the energy and space differential operators (see these operators 
in Eq. H]for 5 = const and Eq. [6] respectively). <p(x) is the occupation number of the seed 
photons, and s(r) their spatial distribution in the bounded medium. The solution of equation 
(jZj) can be conveniently expressed as a series n(x,r) = ^2ckNk(x)Rk(r), where i?fc(r) are the 
eigenf unctions of the space operator L T and Ck are the expansion coefficients over the seed 
photon spatial distribution s(r). We calculate the Comptonization spectrum using only the 
first term k=l of the series (TMK97), which is obtained solving the equation 

, 2 d 2 N(x) . 2 A , . . dN(x) 
fbx 2 „ K + (x 2 + Af h x + 5x) 



d 2 x ' dx 

<f(x 



+(4x - i)N(x) = -7- , 



where 



is given for \i = 2, r e jj =< r > (see Eq. |5]) and 7 = Here is the first eigenvalue of 

the space operator, L T specifically 

L t R(t) + X 2 R(t) = 0. (10) 

For the particular case (V^/c) 2 36, /b ~ 1 and it is possible to express the analytical 
solution of equation (|SJ) in terms of the convolution of the seed photon spectrum (p(x) with 
the Green function (GF) of the energy operator L x according to 



00 



^tb(z) = / G(x,x )(p(x )dx , (11) 

where the GF is defined as 







G ^ X ° ) = x^k)< 



( / \ a+3-<5 
X \ 

iFx(a, q, x)I(xq, a, 5), for x < xq; 



x 



a 

X 



'12) 



— ) iFi(a,g,z )/(x,a,5), for x > x , 
x 
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where 



I(x, a, 5) = e t {x + t) 



a+3— (5 j.q— 1 



(13) 



Jo 



and 



q = 2a + 4 - 5, 



(14) 



while i-Fx(a, b, z) is the Kummer confluent hypergeometric function. Note the change of sign 
of S with respect to the same GF reported in F08 because here we are considering a diverging 
flow. The energy spectral index a of the GF is given by 



and the normalization Cn = a (at — 5 + 3) is chosen in order to conserve the total photon 
number, namely 



It is worth pointing out that oc 1/tq, where r is the total Thomson optical depth 
of the shell, so that 7 <C 1 corresponds to Otq ^> 1, namely for the case of saturated 
Comptonization. In the case of a pure thermal motion (5 = 0) this condition implies a ~ 
(see Eq. |15j). and the emergent spectrum obtained from equation ( TTTj) has a Wien-like 
shape, no matter which is the spectral distribution <p(x) of the seed photons. On the other 
hand, if bulk motion is present, then a ~ 5. 

In Figure [2] we present a typical example of the shape of the GF for different values 
of the bulk parameter 6 defined in equation ([9]), for the case of soft photon monochromatic 
injection (xq = 0.1) and 7 < 1. While in Figure |3] we show a convolution of the Green 
function with a BB seed spectrum (see Eq. [TT]) . It is evident that the important role is 
played by bulk Comptonization in diverging flow in determining the slope of the high-energy 
wing and the position of the high-energy cut-off, both quantities modify the position of the 
peak in the EF(E) spectrum (see Fig. 1 in F08 for comparison with the case of converging 
flow). 

As already outlined above, these analytical results can be obtained for the case /b=l 
(see Eq. |SJ). On the other hand, if /b > 1, the equation has to be solved numerically, and 
we report the details of the numerical solution in Appendix |A] 




(15) 




(16) 



-9 - 



2.3. High energy power-law component in the GRB prompt emission 

The solution of equation (jSJ) clearly shows that above the energy peak E p in the EF(E) 
diagram, the spectrum has an exponential rollover (see Fig. [3]). This spectral feature is 
characteristic of Comptonization both in a pure thermal or subrelativistic moving cases. 
For a thermal Comptonization, the cut-off energy E c is dictated only by the plasma electron 
temperature, with E c >2kT e , while when bulk motion is present (inflow or outflow), the bulk- 
parameter S also concurs in dictating the i^-value (see Figs. [2] and [3]). In both cases, the 
high-energy exponential rollover arises from the presence of the thermal recoil-effect oc N(x) 
in the radiative transfer equation modified by the term oc 5 (see Eq. H]). 

However, it is a well established observational result that most (albeit not all) of 
the observed prompt GRB spectra show an extended PL above the peak energy in the 
EF(E)— diagram [where F(E) is the energy spectrum], which is phenomenologically de- 
scribed by the second index T 2 of the Band function. This high-energy feature cannot be 
reproduced by the solution of Equation (jSJ), and a second component to describe the broad- 
band spectrum is required. One possibility is to simply add a PL to the thermal spectral 
component, as was done, e.g., by RP09, who fitted the BATSE spectra with a BB+PL 
model. From a physical point of view, this model can imply a different origin for the appar- 
ently thermal and non-thermal part of the GRB spectrum. It is also well known that the 
inclusion of a simple PL in broad-band spectral model can introduce some biases, because 
this component, intended to be used for describing only part of the spectrum (usually at 
high energies), actually covers the whole spectral range, producing undesired effects such 
as overestimation of the flux at low energies ( the effect getting strong er for steeper photon 



index T) (Frontera et al. 2012, in preparation; iGhirlanda et al. 1 120071 ) 



Here, we investigate a different possibility; there is a general theorem that the solution 
of an equation with linear differential operator and a source term (such as Eq. [8]) is ob- 
tained by the convolution of the GF of the given operator with the source function itself, 
according to Eq. (11 II) . Examples for this convolution, in addition to the results of previous 
section, are reported for the case of pure thermal Comptonization (ST80), the rmal plus bulk 



Comp tonization (TMK97, F08), and Comptonization in strong-magnetic field (IBecker et al 



20071 ). The resulting energy spectrum can be written in the general form as 



F{E) = j^jHE) + A x <p(Eo) * G(E, Eq)], (17) 

where the two terms on the right-hand side represent the direct and Comptonized (convolved) 
part of the seed spectrum, respectively. The weighting factor A is related to the physical and 
geometrical configuration of the system. The high-energy behavior of the spectrum depends 
on the kernel G(E, Eq); for example, if it is a broken PL with no cut-off, then the convolution 
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with <p(E) (see Eq. HD produces a PL (BMC model in the XSPEC package). 

Thus the observed GRB prompt spectra, phenomenologically described by the Band 
function, can be obtained by an equation of the form (1171) with A> 1, the source term (p(E) 
obtained by the solution of equation (jSJ) and using a broken-PL GF with no recoil-term. The 
latter can be conveniently described by the following form: 

{/ \ a b +3 
/ < -Ob 
If) for x > x , 

where the subscript b has been used to avoid confusion with the index a of the GF related 
to the thermal plus bulk energy operator (Eq. [T2|) . Note that this is the GF of the energy 
operator in equation flU) with fh=l, 5=0 and the recoil term oc x A N in square brackets 



dropped ( jSunyaev &: Titarchuk Ill980l ). which leads to the equation 



^ 2 oHy ^ ^ X ^L - a h (a h + 3)W(x) = -a h (a h + 3)ip(x), (19) 
ax z ax 

where <p(x) = N(x) is obtained by solving equation (JHJ). However, the spectrum obtained 
as a solution of equation flH]) has no analytical representation, so performing a numerical 
convolution according to equation (ITT]) where Gb p i(x, x ) is given by equation ( ITS]) and <p(x) 
by N(x) in equation (jHJ) can be time-consuming in terms of CPU computation. We thus again 
numerically solved equation ( IT9|) using a finite-difference method. From the mathematical 
point of view, there is perfect equivalence between the analytical solutions given by the 
convolution above described and direct numerical solution of the equation. In our model 
is a free parameter, related to the high-energy PL photon index T 2 of the Band function, by 
a b = -T 2 - 1. 

To better explain this behavior, in Fig. H] we show two spectra obtained from the 
BMC model of the xspec package (TMK97), where the source function <p{E) is assumed 
to be a BB occupation number. Note that when A ^> 1, only the second (convolution) 
term contributes to the spectral formation, and the low-energy and high-energy parts of the 
spectrum smoothly join at the peak energy, in a qualitative way which is very similar to the 
Band function. 

This is actually the approach we adopted to develop our model: we first solved the 
Comptonization equation in the presence of subrelativistic bulk motion (see Eq. [8]) and then 
we convolved it with the broken-PL GF of equation (IT8|) . setting A ^> 1 in equation (ITT]) . In 
Appendix [B] we describe the numerical methods used to address these issues. 
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3. The Amati's relation and its interpretation 

In spite of the still open discussion about the impact of possible selection effects on the 
correlation between the peak ener gy of EF(E) diagra m E p and the GRB radiated energy 



or luminosity (Amati relation, AR, lAmati et al. 1120021 ) . it is a matter of fact that all GRBs 



with k nown redshifts, but one (GRB 980425), follow this relation [see, e.g., lAmati et al. 



( 120091 ) and references therein]. In additi on, it is found that also within single GRBs, this 



relation holds (Frontera et al. submitted; iGhirlanda et al. II2010I ). 



Thus it is crucial to understand the origin of this relation, namely which is the mech- 

1 /2 

anism (at first glance universal) which gives rise to the relation E p oc E-J Q where E p and 
E iso are the peak energy of EF(E) diagram and isotropic radiated energy (fluence) of GRB 
respectively during the prompt phase. 

The main issue to be investigated is to check how the parameters of our model concur 
in determining the energy peak value E p in the EF(E) diagram and the total luminosity 
Lg r b which integral over prompt emission time is Ei SQ . 

In Fig. [5] we show the resulting spectra obtained by fixing the parameters of the early 
sub-relativistic phase and changing the spectral index «b of the GF defined by equation (TT%j) . 
which is responsible for the high-energy PL component. Given that E p very slightly depends 
(within a factor less than two) on different ab-values, we fix at, about 1.5. In Fig. |6]we show 
the theoretical EF(E) spectra obtained using our model alternatively changing /3o, kT e , 5 and 
fcTbb, namely the parameters which characterize the subrelativistic phase. In addition we also 
show an example of the spectrum obtained using the solution of Eqs. flE]) and f lT9"j) compared 
them with the Band function in Fig. [7J The result is that the leading quantity in influence of 
E p is the electron temperature kT e . This result is not surprising because it is evident that in 
any photospheric model with spectral emission dominated by Comptonization, the electron 
temperature represents a key ingredient, both in a pure thermal and dynamical cases (for 
|3 < 1). On the other hand, it is worth pointing out that E p is independent of values of the 
BB-like seed photon temperature kT hh . Thus different values of kT hh do not change E p , but 
they determine the total (resulting) luminosity, as we discuss below. 

The question that naturally arises is whether t he observed depend ence of E p on E [so 



(or on luminosity -L gr b) and its intrinsic dispersion (lAmati et al. 1 120081 ) are the result of a 
multi-parametric combination of the subrelativistic phase plus the effect of the non-thermal 
process or is the fundamental effect of 7— ray emission. 

As an example, in Figure we report the results of the fit using our model (see Appendix 
B) of the time averaged prompt emission spectrum of GRB 990705 obtained with BeppoSAX 
and, for comparison, with the Band function. With our model, the best-fit parameters 
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are kT bh = l.8±\% keV, kT e = 94+f° keV, 5 = 0.32±°;^, a b = 1.50lg;$ and photospheric 
radius R ph = 1.5 x 10 13 cm, with x 2 /dof=168/182. We fixed 7 = 5 x 10~ 3 (as a 7-value 
strictly equals to zero produces computer underflow) and outflow velocity (3 = 0.2. The 
Band function on the other hand provides a low energy photon index T 1 = -0.86+ 3 ;^ and 
E = 292±^| keV with x 2 /dof=167/182. In this case however the high-energy index is not 
constrained and we fixed = —2.5 (note that because of sign convention «b = —1^2 — 1)- 

The small values of the best fit parameters of our model 7 and 5 correspond to the 

case of the saturated Comp t onization for which the resul ting index a <C 1 (see Eq. [TBI). 

Sunvaev &: Titarchuk I (jl980l ). ISunyaev fc Titarchuk Ul985h [and see also lChakrabarti fc Titarchuk 



( 119951 )] derive the formula for the Comptonization enhancement factor rj comp which is a ratio 
of the resulting luminosity, that is in our case, the GRB luminosity L grb to the injected 
luminosity of soft photons L soit . Namely 

r/comp = = q X0 (a)x^\ (20) 

-^soft 

where 

= a(a + 3)r( ; (2 + a Tir )r(1 " a) (1 " xl ~ a) ' (21) 

and x = 2.7kT hh /kT e . Thus when a C 1 

^comp oc kT e . (22) 

But the flux of soft photons illuminating the hot spot (Compton cloud) L so f t is 

L S oit = ^-Bsoft'S'hs (23) 

where B so f t is the intensity of the blackbody radiation of the star and Sh s is the surface area 
of Compton hot spot. 

The thermal wave propagates in the hot spot with plasma velocity V p whose value can 
change from one GRB to another. 

As a consequence for each GRB Sh s oc (V p t) 2 and then 

S cc ex V p 2 oc kT p = kT e . (24) 

Thus the luminosity of the GRB hot spot L gr b should be 

-^grb = ^comp-^soft OC (kT e ) 2 . (25) 
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In order to calculate the GRB fluence £^i SO ,grb one should integrate L gr b over the GRB 
prompt emission time T pr , namely 



(26) 



If the time-scale of the GRB prompt emission and its shape is more less the same for any 

1/2 * [ 

/ iso which is precisely seen in the Amati relation [Amati et al. 



(2002 


) and 


i p ^ | 
Amati 


'-■-pi ^ 
(2006)] 



1 12 

correlation E p oc E^ Q can be caused by the spread of the parameter values that characterize 
the outflow evolution of each GRB like the bulk velocity V& and the spectral index of the 
relativistic phase a& (see Eq. [18] and Fig. [7]). Another reason for the spread in the Amati 
relation is the different temperature of the seed (soft) photons, which is related to the 
emission of the star itself. 



4. Temporal variability 



One of the most established observation al results of GRBs s ince the beginning is the 
high temporal variability of their light curves (IMeegan et al. Ill992l ). which can extend down 



to 10~ 3 s. More recently, in addition to the short time variability, also a slow time variabil- 
ity, with time scales from several seconds to ~ 100 s has been discov ered, with the GRB 
prompt light cur ves being the s uperposition of both components (e.g., IVetere et al. 112006 



Gao et al. II2012I ). In particular. IVetere et al. I (120061 ). by performing the systematic analysis 
of a sample of GRBs detected with the WFCs aboard BeppoSAX also showed that the slow 
component is more pronounced at lower energies (2-10 keV) than at high energies (10-26 
keV). Accumulating evidence from the Swift era suggests th at at least the slow component 
of the GRB variability is closely related to the central en gine Liang et al. Ibood ). which was 
confirmed also by other investigations (IGao et al. 1 120121 ). 



Now we demonstrate that, in our scenario, the slow component with temporal variability 
of order 1-100 seconds is related to the Compton cloud while the fast variability component 
with scale of <C 1 s is presumably related to the relativistic outflow (jet). In order to estimate 
the time scale of emergent radiation from the hot Compton cloud one should take a difference 
of the arrival time from the outskirt t out and center of Compton cloud t c (see Fig. [[]). Namely 

At = t out - t c , (27) 

where t out = p ou t/c, t c = d/c and d, p out are the distances from the observer to the center 
and outskirt of the hot spot respectively (see Fig. [p. Thus 



Pout 



R 2 sinV+ [d+ {R-Rcos<p)} 2 



(28) 
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where R is the radius of the star and ip is angle between the vector directed from the center 
of the star to outskirt point and vector directed towards the observer. Then for R <C d we 
find that 

At = t out -t c faR(l- cos ip) /c (29) 

and for (p <C tt/2 we obtain 

At « p 2 R/2c. (30) 

For example if p < 0.1 (or ip < 5° in degrees) then the variability time scale of the 
emergent radiation observed from the Compton cloud should be 

At < 15(.R/10 14 cm) s. (31) 

which is a typical long time scale of the GRB prompt emission (see above). 

On the other hand short time scale variability (t < 1 s) is presumably related to 
relativistic outflow (see section 12.31 and below). Variations of X-ray— 7 radiation in the 
outflow in time scale T detected by the Earth observer as 

T = T[l-{V b /c)cos6) (32) 

where V& is the outflow velocity in the relativistic phase and 9 is angle between the direction 



of the outflow motion and the direction to the Earth observer [see iRybicki &: Lightman 



(119791 )] . If the observer sees the outflow motion at 6 <C 1 then the time scale is very short 

T = T/(2r£) (33) 

where is the Lorentz factor of the outflow(jet). It means that the short variability time 
scales should be only seen when the jet is directed towards the observer. The short temporal 
variability seen in GRB of order of 10~ 3 — 10 _1 s can be considered a signature of the 
outflow for which time scales of X-ray— 7 variations of order of seconds are seen by the 
Earth observer as 10 -3 — 10 _1 -s variations. In other words the long time scales of order 
10-100 s (intrinsic variabilities of X-ray radiation in the hot spot region) are related to the 
Comptonized radiation of the hot spot in the subrelativistic phase while the short time scales 
of order 10~ 3 — 10 _1 s, as seen by the Earth observer, are originated mainly in the re lativistic 



outflow phase similar to that in the standard fireball scenario (e.g. iMeszaros II2006I ). 



Using the emergent spectra itself one cannot distinguish the contribution of Compton 
cloud and the outflow in the spectrum below the peak energy E p because the Comptonization 
spectrum is wider than the upscattering Green function of the outflow, at least for at&>1.5 
(see section 12.31) . The short and long time variabilities should be seen in the same energy 
range at E^E p but the short one should be mostly observed alone at high energies (above 
E p ) where the pure extended power-law component is detected. 
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5. Discussion 

We have developed a Comptonization model aimed to describe the prompt spectral 
emission of GRBs over a broad energy range. The assumption is the presence of two ingredi- 
ents: a thermal bath of soft X-ray photons and the existence of a hot plasma outflow which 
should be fo rmed as a con s equen ce of the explosion of a star. Both ingredients were also 



assumed by lLazzati et al. I (120001 ) for their Compton-drag model to interpret the prompt 
emission from GRBs associated with supernovae. In our model, the hot plasma outflow 
has a bulk velocity which is initially subrelativistic (Vb/c < 0.1 — 0.2) and later it becomes 
relativistic outside the star photospheric radius. Note that a mode l consisting of an early 



subre lativistic followed by a relativistic phase has been proposed by (jLyutikov fc Blandford 



20031 ) in the context of GRBs originated by electromagnetic outflow. 



Our model consists of two parts: the first one describes the spectra up to the peak energy 
in the EF(E) diagram and is obtained from the solution of the Fokker-Planck expansion of 
the radiative transfer equation which takes into account multiple scattering of a BB-like 
seed photon population off a hot and thick electron plasma [Compton cloud (CC)]. The 
relatively soft BB population (kT^t 1 keV) presumably originates in the relatively cold 
star photosphere and illuminates CC (see Fig. [p. The spectral output of this early phase 
is characteristic of saturated Comptonization, but the shape slightly deviates from a pure 
Wien law because of the modification due to outflow motion (see Figs. HJandE]). In this 
case, the peak of the spectrum in the EF(E) diagram falls slightly below E p ~ 4fcT e , with 
respect to that for the standard Comptonization where the bulk outflow velocity Vf, — 0. 

The second part of our model describes the prompt emission spectra beyond the peak 
energy E p . This was obtained by an assumption that, beyond the photosphere, the outflow 
velocity becomes relativistic and provides a further upscattering of the already Comptonized 
spectrum. We discuss now some relevant points related to our model. 

It is important to estimate the thermalization optical depth r e g, which also appears in 
the definition of 5 (see Eq. [5]). This quantity indeed is not an explicit parameter of the 
Comptonization equation (JS}. Using the best-fit values of kT e and 5 for GRB990705 (with 
A,;: 0.2) we find r cS ~ 2. 

From the continuity equation in spherical symmetry we find 

N « 4 x 10 14 ( m//1 ° ^ cm" 3 (M) 

(Wo.5)r? 4 (/yo.i) cm ' m 

where m = M/M Q (M « 10 33 g s" 1 ), r u = i?/(10 14 cm), £ n = 2vr(l - cos 0) and 9 is the 
half opening angle of the subrelativistic outflow. 
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Then optical depth of corona located on the top of the star is 

T ° = ~ 7 (&/0.5)r» M (/Jb/0.1) ' <35) 

where if is the thickness of the hot corona. 



5.1. The resulting GRB luminosity 

The total luminosity of the observable GRB emission can be estimated analytically 
using the same procedure defined in ST80 and CT95. Let us first define the BB luminosity 
illuminating the hot Compton cloud as 

L bb>il = 5 x 10 46 (^/0.5)(A;T bb /0.1 keV) 4 r? 4 erg s" 1 (36) 

Using the formula for the enhancement factor rj comp (see equation [20]) for the case of saturated 
Comptonization (a <C 1) we obtain 



= V comp L bb ,u ~ 5 x 10 49 (T 1 ( t/r (b 1 b) )(en/0.5)(T (b 1 b) ) 4 r 1 2 4 erg s"\ (37) 

where T$ = kT e /(100 keV) and T (bb) = fcT bb /(0.1 keV), respectively 

If the high-energy power law index of the hard component, which is presumably formed 
due to the inverse Compton effect, is a b >2 then the related enhancement factor r] re \ 

_ a h (a h + 3) 

^rel - 7 Try . 7T (38) 

(a b - l)(a b + 4) 

is order of 1. Thus the resulting enhancement factor r] res ~ r] comp and the GRB resulting 

i • • ; T T comp T comp 

luminosity L grb = r] rel L grb p ~ L grb p . 



5.2. On the origin of the high-energy power-law component 



As we already outlined in Section 12.31 subrelativistic bulk motion Comptonization is 
characterized by a rollover feature around an energy which is of the order of the electron 
temperature, slightly modified by the presence of matter velocity field. In any case, this 
proce ss cannot give rise to any PL component extending up to GeV as observed in GRBs 



Abdo et al. I Il992l ). Actually a second component must be included in the X/7-ray 



(e.g- 

spectral model. If the subrelativistic phase Comptonization spectrum and the high-energy 
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PL have different origin, then in terms of spectral modeling the latter could be described 
by simply adding a PL spectral feature. However, a physically unconnected origin of the 
two components seems unlikely, given that all the GRB spectra show a smooth joining of 
the low-energy and high-energy features around the peak energy, and the Ge V emission for 



most Fermi GRBs is consistent with the extrapolation of the MeV component ( iZhang et al. 



2011, hereafter Zll) 



Here, we suggest this physical link in terms of the solution of Equation (lT9"j) where 
the bulk Comptonization spectrum of the subrelativistic phase (Eq. |H]) is convolved with a 
broken-PL GF (Eq. [18]) leading to the emergent broad-band spectrum (see Figs. EH?]). 

The critical point to be addressed is of course the nature of the upscattering process 
that gives rise to the suggested GF (Eq. IT8"j) . The apparently non-thermal origin of the high- 
energy GRB component leads to include mechanisms of particle non-thermal acceleration in 
shock regions. The standard fireball model actually allows the possibility to consider shocks 
as possible sites for particle acceleration, which eventually lead to non-thermal distributions 
N e (E) oc E~ p . In this context, the high-energy PL would be produced through synchrotron 
and/or inverse processes by this non-thermal particle bath. An intermediate case where 
the electron population has a Maxwellian distribution plus a supra-thermal tail was also 



considered by iTavani I ( 119961 ) . 



One of the key assumptions of the fireball scenario is that the GRB outflow is highly 
relativistic since the beginning of the event, with T > 100. In our model, however, the 
first phase of the GRB phenomenon is subrelativistic. How to conciliate this significant 
difference? Is really an unavoidable requirement that the photosphere is ultra-relativistic? 
Let us consider this possible scenario: when the top bounce of the hot shell emerges from 
the star photosphere with bulk velocity V^/c < 0.1, the underlying star photospheric lumi- 
nosity is huge (L* > 10 10 L Edd ), and it illuminates the shell itself, which subtends a solid 
angle Q/Att < 1. This particular configuration (hot plasma, anisotropic high radiation pres- 
sure) could be suitable to give rise to an efficie nt plasma acceleration through the so-called 



Compton- rocket effect" (jCheng fc Odell Ill98ll ). For instance, in the case of an illuminating 



point source, the ratio of the r adiative forc es acting on a hot and cold plasma, respectively, 



is /hot //mid = 2/3 < (7/3) > f lOdell 11198 lh . where < (7/?) > is the quad ratic average mo- 



mentu m of a Maxwellian electron distribution. Numerical simulations by iRenaud fc Henri 



( 120001 ) show that electron-positron jet populations with non-thermal distribution illuminated 



by an underlying accretion disk can be accelerated up to Lorentz factors of several tens. 

It is beyond the scope of our paper to investigate quantitatively the hydrodynamical 
configuration leading to plasma acceleration under the underlying star radiation field. Here 
we just claim a qualitative scenario where the accelerated plasma may interact with the 
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interstellar medium and/or star wind producing shocks with associated particle non-thermal 
populations. The high-energy photons of the subrelativistic phase then interact with these 
electrons via inverse Compton effect (Fig. [T]) and the GF of such a process in given by 
equation (|18p . Note that in this case, the radiative production mechanism is not qualita- 
tively different from the standard fireball model. No matter whether synchrotron or inverse 
Compton dominate the spectral formation, what is important is the total electron Lorentz 
factor, r tot ~ r sh 7 e , where r sh and 7 e are the Lorentz factors of the expanding shell and 
of the electron in the shell frame, respectively. A quantitative treatment of the radiation 
transfer problem is needed in this case. 

Other models have been proposed to explai n the GRB prompt emission in the fram ework 



of inverse Compton processes (see reviews by IMeszaros 1 120061 ; iGehrels et al. II2012I ). The 



key difference between our model and the other photospheric or synchrotron scenarios is 
that in our model case the energy peak of the EF(E) diagram is independent of any Lorentz 
r l— factor associated to relativistic expansion. The main parameter driving the peak position 
(the electron temperature kT e ) is essentially the same in the fl uid and observer f r ame ( apart 
from cosmological redshift effects), unlike that suggested by iGhirlanda et al. I (120121 ) who 
claim that the extension of E p energy of the EF(E)— diagram to a few MeV is mostly due 
to distribution of Lorentz factors in the relativistic shell (or jet). 



5.3. The scenario in the light of the Fermi results 



The thermal nature of GR B spectra is still subject to debate. Blackbody components 



are p r edicted in various m odels ( IMeszaros fc Rees 112000c iPeer et al. Il2006t iThompson et al. 



20071 ; iToma et al. 1 1201 if ) and actually they have been found in time resolv ed GRB spectra 



using the BeppoSXX Gamma- R ay Burst Monitor and Wide Field Camera 2 ( IFrontera et al. 



20011 ) and BATSE instruments (IGhirlanda et al. 1120031 : iRvde 112004 12005L RP09). However, 
from the recent broad band observations with the Fermi GRBM plus LAT instruments it is 
not so clear whether the existence of a thermal component is a general property of GRBs. 
Indeed, on the basis of the systematic and detailed time-integrated and time-resolved analysis 
of a sample of 17 GRB data observed with GBM and LAT, Zll found that most of the GRBs 
in this sample could be fitted with a Band function over the whole Fermi energy range. Only 
in two cases (GRB 090902B and possibly GRB 090510) the spectrum was peculiar, in that it 
could be described by a BB+PL model, similarly to what reported in RP09. The paradigm 
of Zll is that the Band-only spectra are better consistent with a Poynting-dominated rather 
than baryon-dominated flow, thus discarding in most (albeit not all) photospheric 
scenario. Against a pure photospheric emission the authors argue that the extension is up 
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to GeV energies. 

The preponderance of the non-detection of high temperature BB-like features in the 
Fermi GRB spectra as reported in Zll is not, in our opinion, necessarily in contrast with 
our photospheric scenario. Indeed, in our proposed model the seed thermal photons and 
the electrons are decoupled and the low-energy slope of the emerging spectrum (up to the 
peak energy E p ) is mostly dictated by the optical depth and plasma temperature of the hot 
corona and slightly modified by subrelativistic outflow velocity Vb <C c (see Figs. EJandE]). 
The Band function itself can thus be yet the result of photospheric Comptonization process, 
but with a different configuration of the radiation field and electron plasma (decoupled) 
with respect to the standard scenario (where both are coupled). It is also worth noting that 
the above-mentioned BB features detected in BASTE and Fermi spectra have characteristic 
observer-frame temperatures of orders of hundreds keV, and thus must not be confused with 
the cool BB seed photons (<1 keV) of our model. More specifically, we identify these hot 
BB component with the optically thick subrelativistically expanding plasma of the earlier 
phase in our scenario. 

Actually, an important theoretical prediction of our model, to be tested with observa- 
tions, is the spectral steepening at low-energies in comparison with the Band function (see 
Fig. [7]). The change of slope is predicted to occur below ~ 4fcTt,b keV in the EF(E) dia- 
gram, where kT^ is the seed photon BB temperature, and it does represent the low-energy 
tail (Rayleigh- Jeans law) of the input BB. Thus a spectral slope F(E) oc E 2 should be a 
universal GRB observational feature at soft X-ray energies (see also Fig. E]). Unfortunately, 
this issue cannot be solved by Swift or Fermi as their energy thresholds are about 15 keV or 
8 keV, respectively. The BeppoSAX GRBs observed with both GRBM and WFCs are more 
suitable to be used to test our model, even if the lower threshold of the passband (2 keV) is 
nearby the limit of the energy range foreseen for the seed photons. A test of our model has 
been performed using the BeppoSAX data(Frontera et al. 2012, in preparation). 

6. Conclusions 

We have developed a spectral model aimed to describe the broadband prompt emission of 
GRBs. We propose that the spectral emission during the prompt phase, phenomenologically 
modeled by the Band function, is the result of an earlier phase where soft BB-like photons 
are Comptonized by an optically thick and hot electron shell (T e ~ 10 9 K), something like 
a Compton cloud sub-relativistically moving outwards the star surface. On the other hand 
in the relativistic phase, these Comptonized photons are subjected to a second upscattering 
process which can be mathematically described by a broken power-law Green function whose 
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spectral index models the high-energy slope of the Band function. An important prediction 
of our proposed model is that the peak energy in the EF(E) diagram originated in the early 
subrelativistic phase (see Fig. [Q is directly related to the plasma temperature of the hot 
plasma T e . We demonstrate that the resulting luminosity of X/j-rays luminosity of GRB 
Lg r b is proportional to (kT e ) 2 (see Eq. 125]) . 

In fact , Lg r b is a product of the Comptonization enhancement factor 77 comp and lumi- 
nosity of the soft blackbody photons L^u but in the case of the saturated Comptonization, 
when a<l, f] corav , oc kT e (see Eq. [20]) but Lbb,n ( or Araft) is a l so proportional to kT e because 
the surface area of Compton cloud illuminated by the soft photons is proportional to kT e 
(see Eq. I21E]) . Thus we claim that the model dependence of L gr b oc (kT e ) 2 on the hot plasma 
temperature T e explains the observed Amati relation in which E iso = J L grh (t)dt oc E 2 . It is 
worth noting that the peak energy E p of the emergent Wien spectrum should be equal to 
3kT e . 

Our model also explains both the observed long and short temporal variability during 
GRB prompt emission phase. The former one is related to time scale of the radiation formed 
in the sub-relativistically moving Compton cloud while the short time scale is related to time 
scale of radiation formed in the sub-relativistically moving Compton hot spot while the short 
time scale is related to the time scale of radiation coming to the Earth observer from the 
relativistically driven outflow. 

We do not quantitatively investigate the physical conditions underlying the origin of 
the high-energy component. Different scenarios, not necessarily excluding each other, can 
be operating such as a second inverse Compton process off non-thermal electrons accelerated 
by super Eddington radiation pressure (for which L ^> LEdd) as the hot shell moves out of 
the star. Self-synchrotron Compton or direct synchrotron effects can not be excluded. An 
important prediction of our model is that the peak energy in the EF(E) diagram originates 
in the early sub-relativistic phase and is proportional to plasma temperature kT e and the 
resulting luminosity L gr b is proportional to (kT c ) 2 . This dependence is the same, after 
cosmological corrections, in the source and observer frame. In fact, no fine tuning related to 
some Lorentz T-factor of relativistic expansion is required. We claim that this result of our 
model naturally explains the physical origin of the Amati relation. 

A systematic application of our model to a sample of time-resolved GRB prompt spectra 
of the BeppoSAX archive by Frontera et al. (2012, to be submitted) proves that in all the 
cases the x 2 -values are equal or better than those derived with the Band function. To our 
knowledge, this is the first time that a model for the GRB prompt emission other than the 
Band function has been developed under the X-ray spectral fitting package XSPEC. 
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A. Numerical solution of Fokker-Planck radiative transfer equation using 

finite differences 

We consider the Fokker-Planck expansion of the radiative transfer equation for the zero- 
moment occupation number N(x), in the presence of subrelativistic bulk motion (see Eq. |H]): 

„ 2 d 2 N(x) . O „ , f S dN(x) /A \ -AT / \ ¥{ X ) /A-\ 

Kx 2 K ' + (x 2 + 4f h x + 5x)—^- + (Ax - 7 A x = -^V, Al 
d z x dx x 6 

where x = E/kT e , 9 = kT e /m e c 2 , /b = 1 + 3/3 2 /0, while 7 and 5 are defined in Sect. 12.21 
For <p(x) we assume a BB spectrum (see Eq. IBlj) . 

In terms of the zero-moment intensity J(x) = N(x)x 3 the equation becomes 

fx* + ( x 2 _ 2fx + 5x) + (x _ 7 _ 35) J(x) = (A2) 

In order to provide logarithmic binning of the a dimensional energy x, we introduce the new 
variable y such that x = e y , thus obtaining the new equation in terms of this variable 

+ {eV ~ 3/ + 5) ^lr + (e " ~ 7 " 35) J(y) = ~ v[y) - (A3) 

For numerical integration, we divide the domain of y from 7/ min to y max into N equally- 
stepped sizes h = (y max — y m i n )/N, so that equation can be presented, using finite differences, 

as 

Jn+1 2t/ n -|- Jtx—1 , Jn+1 Jn , T / » \ 

Pn ^ h n ^ h r n J„ = -yj n , (A4) 

where 

p a = f, g n =(e v »-3f-5) and r n = (e to + 35 - 7). (A5) 

The two boundary conditions for the intensity J(x) are 

J[0] = 0, J[A] = 0, (A6) 
resembling the physical conditions J(x) = for x — > and x — )> 00. 
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Collecting the terms with the same index in equation flA4j) and using the boundary 
conditions (IA6[) . we may rewrite the equation as 



a n Ji-l + b n J n + C n J n+ i — <7 n , (A7) 

with coefficients a n , 6 n , c n defined as 



a = 0, b = -1, c = 0, ^o = 0, n = 



a n = -^, b n = - (— + — - r n J , c n = + — , (p n = (p(t n ), for n = 1, ...N - 1{A£ 
on = 0, &n = 1, cn = 0, <^n = 0, n = N. 

For n = Equation (IA7I) can be presented as 

Jo = Jo Ji + -Ko, where Jo = — c o/bo, and -K" = 0. (A9) 

Excluding J from equation flA7j) for n = 1, we obtain 

(a x L + 6i)Ji + ciJ 2 = Pi - cLiKo. (A10) 

Because ai Jo + &i 7^ we also obtain that 

Ji = L x J 2 + K x , (All) 
where we have introduced new coefficients 

Li = , and = — ■ (A12) 

aiL + bi aiL + b\ 

Continuing the process of exclusion of J n we obtain in general 

J n = L n J Q+l + K n , (A13) 



where L n = — - — — , and K n = — - — ^-r~- (A 14) 

a n-^n-l + n 0, n L n _i + n 
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Now, using the right boundary condition J[N]=0 and relation (IA13j) it is possible to find 



J N _ X = L N _i Jn + #n-i, (A15) 

and so on until all Jn-2, Jn-3, - Js, J2, Ji, Jo are calculated. This algorithm was first 
introduced by Israel Gelfand is called by the sweep (progonka) method. 

In order to reproduce the unbroken PL Component of the high-energy spectrum, we 
solved Equation (fT9]) for the zero moment intensity, replacing <p(x) with F t b(x) which gives 

jtfJQ _ 2x dJ^)_ _ + = _ Ft ^ x y (A16) 

a z x ax 



As we have already shown in Sect. 12. 3[ this is equivalent to calculation of the convolution 
integral in Equation (11 ip with the GF given by formula (fl8l) and — ^tb(^)- The 

procedure for the numerical solution of equation (1A16h is identical to that used for solving 
Equation (1A1 [) but just changing the expression for the cofficients a n , b n and c n . 



B. Development of a spectral model for XSPEC 

We develop a numerical model for the X-ray spectral fitting package XSPEC. The model 
can be actually divided into two parts, one strictly physical, in terms of the parameter setting, 
and the second one which is more phenomenological. 

Step 1: Comptonization in subrelativistic outflow 

The first part of the model computes the spectrum obtained as a result of Comptoniza- 
tion of a BB-like seed photon population off a hot electron corona which is subrelativistically 
expanding, according to equation (JTJ. For practical purposes, actually we solved equation 
OH]), which is obtained from the variable separation method, with the bulk parameter 5 
defined in Eq. fl9]). 

The input seed spectrum is given by the Planck's law 

^)=N hh {kT c f e(kT J hh)x _ v (Bl) 

where x = E/kT e , while kT^ and kT e are the BB and electron temperatures, respectively. 

The normalization constant A^b = Rg/ 1 Dm pc i where R$ is the BB apparent photospheric 
radius in units of 10 9 cm and D Mpc is the source distance in Mpc, respectively. Note that 
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this is the source term appearing in the right-hand side of equation (JSJ). We numerically 
solved the latter equation using a fast finite-difference method (see Appendix [X]) rather then 
performing a convolution with GF defined in Eq. (II ip. The free-parameters of the model 
at this step are the BB temperature and normalization (kT^ and N^, respectively), the 
electron plasma temperature kT e , the bulk parameter 5, the 7-parameter (which, in fact, 
is the inverse of the Comptonization parameter Y) and the second power of the (constant) 
outflow velocity f3 2 . 

We note that the outflow velocity /3 appears in Equation (|S]) both explicitly in the 
term /b = 1 + (3 2 /(3Q) and implicitly in the definition of the ^-parameter (see Eq. [9]). 
The linear dependence on (3q via the bulk parameter 5 is related to the effect of Fermi 
first-order photon down (up)-scattering and the net result is photon red-shift or blue-shift, 
depending on the sign of V ■ V b , while the quadratic dependence of (3 in / b provides the 
second-order Fermi process and represents an up-scattering contribution (when averaged 
over angles) term in addition to that coming from the electron thermal velocity component. 
By setting (3q <C 39, it is possible to find the analytical solution for the emerging spectrum 
F th (x) using convolution (ITTj) and the analytical presentation of the GF (see Eq. [T2]l . A 
comparison between the numerical and analytical solution in this limit has shown an excellent 
agreement between the two cases. 

Step 2: high- energy powerlaw component 

We should point out that in principle a general broken-PL GF is described by two 
indexes for E < Eq and E > Eq, while in the GF of equation ffl8l) they both depend on 
Qb- However, for our purposes, it is important that the red wing of the GF, oc E ab+3 , is 
significantly steeper than the blue wing oc E~ ah . Actually, it is this up-scattering term of the 
GF which dominates the convolution process. We show an example of spectrum obtained 
from the solution of Eqs. (jSJ) and (fT9~j) compared with the Band function in Figure [7J It is 
evident their similarity, but we also draw the attention of the reader to the change in the slope 
at low energy of our model with respect to the Band function. The low-energy steepening 
occurs at E < 4 kT^ in a EF(E) diagram, where fcTbb is the temperature of the seed BB 
photons, and it is related to the Rayleigh- Jeans part of the Planckian distribution, where 
F(E) oc E 2 . For the case reported in Fig. Owe assumed fcTbb = 1 keV, so that the steepening 
occurs below 4 keV. We subsequently implemented our model in the xspec package and 
tested it on GRB990705 (z=0.842; Amati 2006), observed by the BeppoSAX/WFC+GRBM 
instruments. 
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Fig. 1. — Schematic view of the explosion of a massive star and associated emergent spectral 
components. At the first stage, the spectrum is formed by Comptonization of blackbody- 
like seed photons (red arrows) placed at the bottom of a hot sub-relativistically outflow 
(SROF or Compton cloud) on the top of the star photosphere (orange region). The emergent 
Comptonization spectrum (bright-blue arrows) is obtained using the solution of Eq. flH]). 
The Compton cloud can be accelerated by the underlying radiation pressure of the star 
photosphere leading to the relativistic outflow (ROF). Inverse Compton scattering of the 
Comptonized photons of the subrelativistic phase off non-thermal relativistic electrons of 
ROF can be the origin of the extended power-law observed above the peak in the EF(E) 
diagram (dark-blue arrows). 
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Fig. 2. — Evolution of the Green function (fl2l) for the specific case of saturated Comp- 
tonization (7 ~ 0) and different values of the bulk parameter 5. Black, red, green and blue 
correspond to 5 =0.5, 0.7, 0.8 and 0.95 respectively. The horizontal scale is defined in units 
x = E/kT e . 
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Fig. 3. — Saturated Comptonization spectra (7 ~ 0) modified by the presence of subrel- 
ativitic bulk motion obtained from the numerical convolution of the Green function (fl2|) 
with a blackbody spectrum of temperature kT^ = 1 keV. Different colors represent different 
values of the bulk parameter 5 according to Fig. [2J A plasma temperature of kT e = 100 keV 
is assumed. 
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Fig. 4. — Emergent spectra obtained as the sum of a BB spectrum plus its convolution 
with a broken power-law Green function (see Eq. [T8|) and with weighting factors A » 1 
(dotted- dashed line) and A=0.1 (dashed line), according to equation (ITTj) . In both cases the 
BB temperature is 1 keV and the Green function spectral index ctb=2. 
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Fig. 5. — Emergent spectra obtained by the solution of Equations (JS} and ( II 9 1) for differ- 
ent values of the up-scattering Green function index «b determined by Eq. ( fl8l) . Other 
parameters of the model are kT^ = 1 keV, kT e = 100 keV, 5 = 0.5, f3 = 0.2, and 7 ~ 0. 
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Fig. 6. — Simulated EF(E) spectra obtained as solutions of Eqs. (JS]) and (fT9|) for different 
sets of the Compton parameters values. In the top left of each panel, the fixed parameters are 
reported. The low-energy steepening of the spectra is the transition towards the Rayleigh- 
Jeans region of the seed BB spectrum. Panel a): (3 =0 (black), 0.1 (red) , 0.2 (green), 
0.3 (orange), 0.4 (blue), 0.5 (pink). Panel b): kT e =50 keV(red), 100 keV (green), 150 keV 
(orange), 200 keV (blue), 250 keV (pink). Panel c): 5= 0.8 (red), 0.7 (green), 0.6 (orange), 
0.5 (blue), 0.4 (pink). Panel d): /cTbb=0.5 keV(red), 1 keV (green), 1.5 keV (orange), 2 keV 
(blue), 2.55 keV (pink). 
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Fig. 7. — Comparison of the EF(E) spectra in the 1 keV - 10 MeV energy range of our 
Comptonization model {dashed line) obtained from the solution of equations ([8]) and f|T9|) and 
the Band function {dotted- dashed line). The parameters of the simulation are: T\ = —1.5, 
r 2 = -2.5 and E Q =350 keV for the Band function, kT hh =l keV, A;T C =100 keV, 5=0.5, 
7 = 5 x 10~ 3 , po = 0.4 and afe=1.5, for the Comptonization spectrum. 
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Fig. 8. — Unabsorbed EF(E) BeppoSAX/WFC+GRBM spectra of the prompt emission of 
GRB990705 (z=0.842; Amati 2006) and residuals between the data and the model in units 
of a using our photospheric model (left panel) and the Band function (right panel). Best-fit 
parameters for both cases are reported in Section [3l 



